Adaptive homodyne phase discrimination and qubit measurement 
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Fast and accurate measurement is a highly desirable, if not vital, feature of quantum computing 
architectures. In this work we investigate the usefulness of adaptive measurements in improving the 
speed and accuracy of qubit measurement. We examine a particular class of quantum computing 
architectures, ones based on qubits coupled to well controlled harmonic oscillator modes (reminis- 
cent of cavity-QED), where adaptive schemes for measurement are particularly appropriate. In such 
architectures, qubit measurement is equivalent to phase discrimination for a mode of the electro- 
magnetic field, and we examine adaptive techniques for doing this. In the final section we present a 
concrete example of applying adaptive measurement to the particularly well-developed circuit-QED 
architecture. 

PACS numbers: 03.67.-a, 03.67.Hk, 32.80.Qk 



I. INTRODUCTION 

The exquisite control demanded by quantum informa- 
tion processing tasks over the initialization, evolution, 
and measurement of physical systems has motivated an 
intense level of study into these processes during the past 
decade. The numerous and diverse proposals for imple- 
menting quantum computing mean that we are learning 
how to control these processes in an array of physical 
systems of varying scale and complexity. 

Several promising avenues for solid state realizations 
of quantum computers have been proposed (e.g. P, 0, 
H, [j, H, 0, A recent trend in such sohd state re- 

alizations, particularly ones based on superconducting 
circuits, has been to couple the qubits that form the 
quantum computer to well controlled harmonic modes 
of the electroma gne tic (EM) field (e.g. circuit-QED Q, 
and Refs. [1, [ol. llOl|V Such qubits-plus-harmonic-mode 
systems closely resemble the cavity-QED paradigm pi- 
oneered in quantum optics (Tl| . with atoms now being 
replaced by the qubits (which can be thought of artificial 
atoms). Typically a system of this type is well described 
by the multi-atom Jaynes-Cummings (JC) model: 
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where is the energy splitting of the jth qubit, ojfis the 
frequency of the EM mode, gj is the coupling strength 
between qubit j and the mode, a {a)) is the annihilation 
(creation) operator for the field mode, and {(jL) is the 
raising (lowering) operator for the jth qubit. We have 
set hto \ and will do so for the remainder of the paper. 
This Hamiltonian accounts for the coherent interaction 
between the elements; additional terms must be added to 
describe dissipation (spontaneous emission for the qubits 



and damping for the field mode) , and driving of the field 
mode (e.g. see Ref. Q). 

In the context of quantum computing, the cavity in 
such solid state cavity-QED systems has a dual function: 
as a bus that conveys quantum information between the 
qubits, and as an enabler of qubit measurement. The lat- 
ter role comes from the fact that the well defined coupling 
between the qubit/s and the EM mode allows one to in- 
fer the state of the qubit/s from the state of the mode. 
In particular, when all the qubits are detuned from the 
cavity frequency, i.e. = \uj^^ — \ ^ Vj, then 
a measurement of the cavity is equivalent to a quantum 
non-demolition (QND) measurement of the state of the 
qubit/s. To see this, we note that in this detuned (off- 
resonant) limit, the dynamics of the system is well ap- 
proximated by the dispersive Jaynes-Cummings Hamil- 
tonian 1121: 
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where Xj = 9]l^j- The usual dispersive JC Hamilto- 
nian has an additional term proportional to o'!^(7^ + 

ctL<5'+ Vi,j i 7^ j that describes the intcrqubit coupling 
enabled by virtual excitations of the cavity. This term 
is reduced in magnitude when the qubits are made off- 
resonant with each other [l3l | , and since this is the regime 
that is most useful for measurement, we will ignore this 
mutual interaction term in the dispersive Hamiltonian. 
Note that the first term in Eq. ^ describes a "pull" of 
cavity mode frequency by the qubit/s, and that this pull 
is dependent on the state of the qubit/s (in the basis). 
Intuitively then, if all the Xj ^^re different (i.e. all the w;^ 
are different), there is a unique cavity pull for each state 
of the qubit/s (in fact, it is not sufficient that all the 
pulls simply be different for a unique determination; see 
section IIII B 21 for the complete condition) . This pull is 
exactly what is used to measure the state of the qubit /s 
using the cavity by virtue of the fact that a measure- 
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mcnt of the phase of the cavity output field determines 
this frequency pull @, For example, when we have 
one qubit, the dispersive JC Hamiltonian takes the form 
H = {ujf + x^z) a^a+ ^{iUa + x)^z, and a probe coherent 
state of the cavity undergoes the following input-ouput 
map: 



if (a,) = 1 



if (a.) 



-1 



(3) 



where 



= tan 



1 is the magnitude of the phase 

change (here, k is the decay rate of the cavity). 

Therefore, in such qubit/s-plus-oscillator systems a 
measurement of the state of the qubit/s is often a phase 
measurement of the oscillator, and the task of qubit/s 
state determination becomes the same as that of phase 
discrimination for a harmonic oscillator mode. Thus, if 
we are interested in improving qubit measurement, we 
are naturally led to the question: how best to realize 
phase discrimination of an EM mode? Quantum state 
discrimination, especially in the optical context, is a well 
studied problem with a long history mostly motivated by 
communication problems |15j . The standard method for 
phase measurement of a field mode is to perform a dyne 
measurement of the mode (l6l . [l7j . This involves mixing 
the mode with a local oscillator on a beam-splitter and 
performing photodetection on the two output ports of the 
beamsplitter; see Figure [TJ A homodyne measurement 
uses a fixed local oscillator phase and thus measures one 
particular quadrature of the input mode, while a hetero- 
dyne measurement uses a rapidly varying local oscillator 
phase, which results in a measurement that samples all 
quadratures of the input mode equally. Both these dyne 
measurements are static schemes - they do not use in- 
formation gained about the unknown phase to update 
the local oscillator phase in order to optimize further in- 
formation gain. Wiseman [31 has shown that adaptive 
schemes can outperform such static schemes for the task 
of phase estimation for field modes. This was followed 
by a study by Bargatin [l^ showing that general coher- 
ent state discrimination is improved by adaptive schemes. 
While Refs. and use different measures of perfor- 
mance, both clearly show an improvement with adaptive 
strategies. Also, Dolinar showed already in 1973 [131 that 
the discrimination of coherent states by photon counting 
can be done more efficiently by using adaptive measure- 
ments and this was recently confirmed with both simu- 
lation and experiment by Geremia et. al. [2l|, [2^. In 
fact, the adaptive strategy detailed in achieves the 
optimally discriminating measurement in the limit of in- 
finitely fast feedback. These observations motivate the 
central question of this paper: can one distinguish opti- 
cal phases, and ultimately qubit states in the cavity-QED 
context, more efficiently usin g a daptive dyne measure- 
ments? The results of Refs. (isl . [lOl [20j would suggest 
so, and we explore the question in detail in the following 
sections. 

Firstly, in section [IT] we provide some background to 
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FIG. 1; (Color online) The typical method for measuring the 
phase of a field mode: a dyne measurement. The input mode 
(for the measurement) is the output of the cavity, and in the 
cavity QED scenario we are considering, it will have an un- 
known phase (0) due to the pull of qubit/s in the cavity. The 
local oscillator (LO) phase, $(t), is constant for a homodyne 
measurement and rapidly oscillating for a heterodyne mea- 
surement. 



adaptive dyne measurements. Then, in section IIIII we 
examine the task of adaptive phase discrimination, first 
looking at the case of two phases (i.e. measuring a sin- 
gle qubit in the cavity QED quantum computing sce- 
nario) and then at adaptive phase discrimination for an 
arbitrary number of phases. In that section, we treat 
in detail an example of adaptive state discrimination in 
circuit-QED [^, a canonical solid-state implementation 
of the qubits-plus-harmonic-modc system. We conclude 
in section HVl with a summary of the results. 



II. ADAPTIVE DYNE MEASUREMENTS 

A typical (balanced) dyne measurement is illustrated 
in Fig. [T] The field mode with the unknown phase {(j)) is 
mixed with a mode in a coherent state of known phase 
(<i>(i)) and large amplitude (/3), and both outputs are 
measured by photodetectors. The output signal, I{t), 
also known as the photocurrcnt, in the interval \t, t + dt) 
is defined as [23, M : 



lit) = 



lim lim 



SN2{t) ~ SNijt) 



= e-*/2(^g-.*(t)+ate'*(*))+^(i), 



(4) 

where 6Ni{t) is the photocount at the ith detector, a,a'^ 
arc the annihilation and creation operators for the input 
mode (with unknown phase), and ^(t) is delta-correlated 
Gaussian white noise that captures the unavoidable shot 
noise of the detectors [l3] • The decaying exponential en- 
velope in the above average signal refiects the fact that 
the probe beam is a finite pulse (we do not consider con- 
tinuous probe beams). Note that by varying $(t) we 
can measure any quadrature of the input mode - for ex- 
ample, <f>(t) = measures the x oc a + a' quadrature. 



3 



PD 1 



cavity 



f ■ N 1"^'*) Lbs 




Phase modulator 





signal 


i 


processing 





P 

Po 








^P2 


Pa 








^^3 





LO 



FIG. 2: (Color online) An adaptive scheme for measuring 
the phase of a field mode. The local oscillator phase is dy- 
namically modified according to the measured values of the 
photocurrent. 



As mentioned above, a rapidly varying samples all 
quadratures of the input mode equally. In the case that 
the input mode is a coherent state, |ae"^), where a is 
real, then the output photocurrent simplifies to 



I{t) ^ 2ae"*/2 cos[$(t) - 0] + ^(t). 



(5) 



An equivalent quantity that will prove to be more useful 
in the following is the photocurrent increment 

I{t)dt ^ 2ae"*/2 cos[$(t) - <\)\dt + dW{t), (6) 

where dW{i) are Wiener increments of mean zero and 
variance dt [13] . Note that the probe signal magnitude a 
functions as an effective signal-to-noise (SNR) parameter. 

Now, consider the modification to the standard dyne 
measurement shown in Fig. [51 Here, the information con- 
tained in the output photocurrent is used to dynamically 
modify the phase of the local oscillator. This introduces a 
feedback loop which opens up the possibility of adapting 
the dyne measurement according based on information 
received in earlier stages of the measurement. 

Wiseman first considered such an adaptive scheme for 
the task of phase estimation. He showed that using an 
adaptive homodyne scheme one can reduce the variance 
of the estimate faster that with any static scheme p^ . 
His result can be understood simply by noting that the 
sensitivity of a dyne measurement, defined as |(5/(t)/(50|, 
is greatest when | sin[(^ — maximized. This can 

be achieved by a such that — = ±7r/2; that 
is, by a local oscillator that is in quadrature with the un- 
known phase. Wiseman's scheme uses this fact to choose 
the local oscillator phase = -I- 7r/2 where 
is the current estimate of the phase (/>. This estimate 
is formed from the measurement current up till time t. 
Thus as the phase estimate improves, the sensitivity of 
the measurement increases. 

We want to use this idea of adaptive measurement to 
discriminate between a given number of optical phases. 



FIG. 3: A phase space diagram showing the possible phases 
in a phase discrimination task with 6 hypotheses. 



III. ADAPTIVE PHASE DISCRIMINATION 

The analysis of phase discrimination can be made 
quantitative by treating it as a formal decision process, 
and more particularly, as a hypothesis testing task. Each 
possible phase forms a particular hypothesis and the task 
is to determine the correct phase from the experimen- 
tal observations: the photocurrent up till time t, which 
we denote by I[o.t)- More precisely, the hypotheses arc 



p = Pj 



one for each 



possible phase, (pj, resulting from a different pull of the 
probe beam frequency by the qubits in the cavity. We 
restrict ourselves to probe beams that are coherent states 
and linear cavity dynamics (so the output coherent state 
amplitude is the same for all qubit states - only the phase 
changes) since this is experimentally most relevant. We 
will comment on non-linear cavity dynamics in the fi- 
nal section of the paper. Figure [3] shows the coherent 
state "constellation" associated with a phase discrimi- 
nation task; the members of the constellation form the 
hypotheses. 

By the likelihood principle [25| , all the information rel- 
evant to the decision is contained in the likelihood func- 
tions for the hypotheses, defined as: 



^3 = Pr(I[o, 



|Wj,'i>[o,t))=Pi-(I[o,t)IP = Pj.*[o,t))- (7) 



As the right-hand side indicates, the likelihood function 
is a conditional probability for the observed measurement 
record given the particular hypothesis. We also condition 
on the local oscillator phase history in anticipation that 
this will become a dynamical variable. The likelihood 
principle asserts that all decisions regarding the hypothe- 
ses are to be made from these likelihood functions and 
therefore we are led to consider them in more detail [l^ . 

Using Bayes' rule, we express the a posteriori proba- 
bilities of various hypotheses in terms of the likelihood 
functions as: 



T'j =Pr(p = p,|I[o,t),$[o,t)) 



Pr(I[o,t)|p = Pj.'^mKj 



Pi-(I[o,t)|$[o,t)) 



AT 



(8) 
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where Q is the prior probabiUty that p = pj, and 
the denominator is a normaUzing quantity: M = 

Pr(V)|$[o.i)) -i:,Pr(I[o,t)|p = A,$[o,t))G =E«AO. 
The prior probabilities are assumed to be known (they 
are the prior probabihties of the qubit/s states), and 
from here onwards we will assume that all qubit states 
arc equally likely and therefore use a uniform prior. In 
this case, note that the a posteriori probabilities for the 
hypotheses {Pj) are simply normalized versions of the 
likelihood functions (Cj). Now, quantum measurement 
theory tells us that 

Pr(%,)|p = p,,$[o,t)) = tr(p,Ft), (9) 

where Ft is the POVM (parameterized by the continuous 
variable t) that corresponds to the observed photocur- 
rent, I[o,t), and local oscillator history, *&[n *v Wiseman 
has evaluated this POVM as being [ll, [11] :' 

Ft = F{Rt,St) = Po{Rt,St)GtiRuSt) 

- Po(i?t,5t)cxpQ5tat2+i?,atj 

X cxp{~a^at) cxp ^^SfO.^ + RfO^ , 

(10) 



where Rt and 5*4 are functionals of I[o,t) that are sufficient 
statistics [25l | which completely capture the influence of 
the measurement record. They are explicitly: 

Rt = i?4I[o,t)] = f e'^<^'^e''/'lis)ds 
Jo 

St = St[Iio.t)]^- I e'2*(^)e-^ds. (11) 

JO 

Note that the second quantity is also a functional of the 
measurement record if the local oscillator phase, $(t), 
depends on the measurement record up till time t, I[o,(), 
as will be the case in our adaptive scheme. PQ{Rt,St) 
in Eq. ()10|) is a normalizing factor which is inconsequen- 
tial to us because it is present in both the numerator 
and denominator of the expression for the a posteriori 
probabilities for the hypotheses, Eq. ([8]), which are ulti- 
mately what we are interested in. Thus we redefine the 
likelihood functions as: 

£j = (ae*"^^! Gt lae*"^^) , (12) 

and evaluate the expectation value to obtain: 



Cj{t) = cxp{-\a\^{l-e-*)}cxp{Re{Staf) + 2Re{Rta*)} 

= Ctexp(-2|a| / (|a|e-^ cos2($(s) - (^j) - e'^/^ cos($(s) - v?j)/(s)) ds }> , (13) 



where Ct = exp{— |ap(l — e *)} is a prefactor that is common to all j. 

I 



Eq. (fT3|) gives us explicit expressions for the likelihood 
functions. Wc will utilize these to form the decision rules 
for the hypothesis testing problem in the following sub- 
sections. Note that the likelihood functions are them- 
selves random because they are functions of the random 
observation process. 

A. Adaptive phase discrimination between N — 2 
phases 

Consider the simplest case of discriminating between 
two possible phases. In the cavity-plus-qubits experimen- 
tal setup this corresponds to having one qubit coupled to 
the cavity mode. In such a case we have two hypothe- 



ses: H+ : p = p+ = jae"^) (ae"^] and TL- : p = p^ = 
\ae -^^) {ae-^^\, for some (p. In a two-hypothesis test, 
the decision policy is simple and obvious (at least once 
the likelihood principle has been adopted). It is based 
on the likelihood ratio, A = which is simply the 

ratio between the likelihood functions for the two hy- 
potheses. The decision policy is to accept Ti+ if A > 1 
and if A < 1 (in the degenerate case of A — 1, a 
random decision can be made) [s^l ■ 

Using the expressions for the likelihood functions de- 
rived above, and assuming equal priors (i.e. both states 
of the qubit are equally likely) we write the likelihood 
ratio as: 
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A(t) = cxp<{-2|ary e-^[cos^($(s)-0)-cos^($(s) + 0)] ds 

+2\a\ I e-"/2[cos($(s) - 0) - cos($(s) + 0)]/(s) ds 



Because of the exponential in this expression, we will 
work with the log of the likelihood ratio; this merely 
changes the decision policy by changing the threshold 
(i.e., InA(t) > or InA(i) < 0). 

Given this decision policy, we can attempt to optimize 
the probability of a correct decision by using an adaptive 
scheme. That is, we want to choose at time t such 
that the probability of a correct decision is maximized. 
This can depend on the photocurrcnt up till time t: 
I[o^t). The probability we want to maximize is: 

Pc{t) = i[Pr(lnA(t) >0 I p = p+) 

+ Pr(lnA(t) < I p = p-)], (15) 



(14) 

I 

where again, we have assumed equal prior probabilities 
for the two phases. This optimization maximizes the 
probability of a correct decision at all times t. An al- 
ternative would be to maximize this probability only at 
some fixed time T, however, we shall find the former eas- 
ier to do and of course, it implies the latter optimum. 



Using the expression for the measurement current, 
Eq. ([6]), we can write this success probability as: 



2p^{t) = Pr(^^ e""/2[cos($(s)-0)-cos($(s) + 0)] dVF(s) >-|a|^ e"" ds[cos($(s) - ?i) - cos($(s) + 0)]' 
+Pr^^ e-''/2[cos($(s) -0) -cos($(s) + 0)] dW{s)<\a\j^ e"" ds[cos($(s) - 



- cos($(s) + 4>)]' 



(16) 



Note that the left hand sides of both probability ar- 
guments are the same; namely, some random variable 
that is a sum of weighted dW{s) terms. The sum is 
not necessarily Gaussian because the $(s) could depend 
on dW{s'),s' < s. However, the left-hand-side can sim- 
ply be treated as a random variable, which we will la- 
bel X(i). Hence the success probability takes the form 
p,(t) = l/2[Pr(X(i) > -xW) +Pr(X(0 < x(<))], where 
X = \a\ /p (is[cos($(s) - 0) - cos($(s) -|- 0)]^. Re- 
gardless of the distribution of X(i) (as long as the do- 
main of the random variable is the entire real number 
line, which it is in this case) this probability will be 
maximized when x{t) is maximized. Furthermore, max- 
imizing xit) is trivial since it is an integral of a positive 
semidefinite quantity: x(t) is maximized when the inte- 
grand is maximized for all s. And it is easy to show that 
e~^[cos($(s) — 0) — cos($(s) + 0)]^ is maximized by the 
choice $(s) = 7r/2 for all s. 

Hence we see that the optimal strategy for distinguish- 
ing two phases is not an adaptive strategy at all, but 
rather a static one that keeps the local oscillator phase 
fixed at 7r/2. The above calculation can be generalized 



easily to the case of distinguishing two arbitrary phases 
ipo and ipi that do not necessarily add to zero: the op- 
timal strategy then is <i>(t) = tt/2 + {ipo + ipi)/2. The 
uselessness of an adaptive strategy in this case can be 
understood by looking at the optimal dyne measurement 
on the phase space diagram of Fig. |4l The optimal 
measurement is such that the local oscillator phase is 
in quadrature with the bisector of the phases that are 
to be distinguished - we will refer to this measurement 
as one that is symmetrically in quadrature. Such a mea- 
surement gathers the most information relevant to the 
discrimination task per unit time. And since the phase 
space is two dimensional, this measurement exists for any 
two phases (po and ipi . This agrees with the recent result 
in Ref. [l^] that the optimal Gaussian strategy (of which 
dyne measurements are a subset) for discriminating two 
coherent state phases is a static homodyne measurement. 

However, if the task is to distinguish between more 
than two phases, such a symmetrically in quadrature 
measurement no longer exists. Therefore in such cases, 
we might expect an adaptive strategy to outperform a 
static one. We will investigate this in the following sub- 
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FIG. 4: A phase space diagram showing the possible phases in 
a phase discrimination task with two hypotheses. As shown in 
the main body, the optimal local oscillator phase is static and 
in quadrature with the bisector of the two candidate phases. 



section. 



B. Adaptive phase discrimination between N > 2 
phases 

When the number of phases to discriminate is greater 
than two, the problem is one of multiple hypothesis test- 
ing. This is a much harder problem and we do not have 
results of optimality here. A sensible decision policy for 
multiple hypothesis testing, when the likelihood princi- 
ple has been adopted and the losses are uniform, is the 
maximum likelihood decision - i.e., choose the hypoth- 
esis with the maximum likelihood function [25|. Given 
this policy, one might formulate a performance function 
for the decision problem as: 



^corr (^) 



i^corr 



(17) 



where Ccorr is the likelihood of the correct phase. Then 
the task would be maximize this performance function by 
a choice of $[0,*)- However, note that the likelihood func- 
tions are non-differentiable functions and therefore many 
of the standard tools of optimization, such the calculus of 
variations or dynamic programming, become unsuitable 
for use. 

Given this situation, we will use the results of the 
last subsection to motivate the following ad-hoc adap- 
tive strategy for modulating the local oscillator phase: 



2 2 



(18) 



where (pM{t) and ipm{t) are the phases with the largest 
and second-largest likelihood function at time t. This 
strategy involves keeping track of the likelihood functions 
for all the candidate phases (all hypotheses) and choosing 
to perform the symmetrically in quadrature measurement 
of the two hypotheses of largest likelihood. There are no 
claims of optimality for this strategy, but it is a sensible 



one to adopt given the decision rule. In the following, 
we will numerically evaluate this adaptive scheme in two 
illustrative situations. 



1. Example 1: a circuit- QED example 

The first numerical study is motivated by the circuit- 
QED architecture for quantum computation d, . This 
architecture is a direct analogy of the cavity-QED se- 
tups in quantum optics and therefore the discussion in 
the Introduction applies exactly. We consider the situa- 
tion where there are two qubits coupled to the strip-line 
microwave resonator that functions as the cavity [28j . 
We assume these qubits have different energy splittings 
and therefore will pull the cavity frequency by different 
amount when the cavity is detuned from both of them. 
Therefore, a measurement of the phase of the cavity out- 
put mode could yield four different values, each of which 
forms a hypothesis: 
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where 



(19) 



tan -^{gf / KAi), and gi and are, respec- 



tively, qubit i's coupling strength to, and detuning from, 
the cavity mode (k is the cavity damping rate, which is 
a function of the capacitive coupling of the stripline to 
the output port). A wide range of values for gi are ex- 
perimentally accessible in the circuit-QED experiments. 
Using the cavity pull reported in Ref. {g'^/nA w 2.5) 
as a guide, we set 0i = 47r/10, 02 = Stt/IO. 

We now compare the adaptive strategy described by 
Eq. to the static strategy of heterodyne detection. 
Static heterodyne detection, where $ is cycled rapidly, 
is known to be superior to the alternative static strat- 
egy of homodyne detection when the number of phases 
to discern is > 2. The average performance of a static 
heterodyne scheme is linked to how fast the phase is cy- 
cled: it improves as the phase is cycled faster but the 
improvement plateaus after a threshold is passed. We 
determined this threshold and operated the static het- 
erodyne scheme above it (numerically, at IOOtt rads/sec). 
We then simulate the measurement current generated by 
a given phase and track the likelihood function of the 
correct phase under the two different strategies. This is 
done 500 times in order to capture the average behavior 
and the results are shown in Fig. [5l 

These figures show the average development of the a 
posteriori probability for the correct hypothesis, Vcorr, 
as a function of time and signal-to-noise, SNR (recall 
that this is measured by a) . Vcorr is of course one of the 



Vj defined in Eq. 



and for the simulation we gener- 



ated a measurement signal consistent with the smallest 
phase, iy9i = 01 — 02 = 7r/10. Hence, Vcorr = Vi in this 
case. The results indicate that, on average, the adaptive 
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(a) Time evolution of the Vcorr for an example 
fixed probe signal magnitude (i.e. fixed SNR, 
a). 
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(b) The Vcorr as a function of probe signal 
magnitude (SNR) a at two fixed times. The 
lines are interpolations between data points at 
integer |o| values to aid the eye. See main text 
for comments on the error bars. 
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(c) The time taken for the Vcorr to reach 0.5 as 
a function of probe signal magnitude (SNR) a. 
The lines are interpolations between data 
points at integer |q| to aid the eye. 

FIG. 5: (Color online) Figures contrasting the evolution of 
the average a posteriori probability of the correct hypothesis 
{Vcorr) under the adaptive and static schemes described in 
section UlI B II In all the figures, the solid (red) line shows the 
evolution under the adaptive scheme and the dashed (blue) 
line shows evolution under the static heterodyne scheme. 

strategy outperforms the static strategy for intermedi- 
ate times. For very short times, the static heterodyne 
scheme performs shghtly better or equally well, because 



the adaptive scheme is ineffective at short times while 
the Bayesian filter accumulates initial data. And at long 
times, the performance of the two schemes becomes sim- 
ilar, although the adaptive strategy converges to large 
likelihood values faster. However, at intermediate times, 
figure [5] shows that the adaptive strategy performs better 
in that: (i) it requires less time for the correct a poste- 
riori probability to reach any threshold value, and (ii) 
it requires a smaller signal-to-noise ratio to achieve the 
same level of certainty. Similar conclusions result from 
comparing the adaptive strategy with a static homodyne 
detection scheme, with the differences being slightly more 
pronounced. It should be noted that these graphs show 
average results; in contrast, a single run of the measure- 
ment would result in an evolution of the observer's state 
of knowledge as illustrated in Fig. [51 Note that in Fig. 
[5jb) the error bars at each data point show the standard 
error of the mean values plotted, i.e., aj^/n^ where a is 
the standard deviation and n is the sample size (these 
error bars are omitted in Fig. Efa) in the interest of clar- 
ity). The small errors bars are assurance that the sample 
size of 500 is sufBcient to capture the average behavior 
well. 




time (sees) 

FIG. 6: (Color online) Example of a single run of the mea- 
surement in the circuit QED example of section UlI B II The 
red (solid) line is the evolution of the a posteriori probability 
of the correct state {Vcorr) under the adaptive scheme and the 
blue (dotted) line is its evolution under the static heterodyne 
scheme. 

Figure[S]shows the results for a particular choice of cor- 
rect state. The average behavior of the adaptive scheme 
when averaged over each of the four possible phase states 
taken as the correct one, and the corresponding average 
probabilities from static heterodyne detection arc illus- 
trated with several examples in Table HI The table shows 
the average conditional probability of the correct state 
for |q;| — 5 and for two choices of time (note that there 
are two averages taken here, one over many runs of the 
task - which results in Vcorr ~ and then another over 
the four possibilities for correct phase state). The table 
clearly shows the advantage of the adaptive scheme, with 
increased probabilities of success at short and long times. 
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Time (sees) 


Detection scheme 


Average Vcon- 


a 


0.2 


Static 


0.8595 


0.0164 


0.2 


Adaptive 


0.9089 


0.0144 


1 


Static 


0.9680 


0.0078 


1 


Adaptive 


0.9951 


0.0026 



TABLE I: A'' = 4 detection: values of the average a posteriori 
probability of the correct phase, Vcorr, and standard deviation 
of this quantity, a, at two time values, one short and one long, 
for fixed SNR (|a| = 5). 



2. Example 2: large N hypothesis testing 

As the number of phases to be discriminated, TV, gets 
larger, the discrimination task approaches the estima- 
tion task considered by Wiseman in Ref. [31 . How- 
ever as the number of qubits in the cavity increases, 
the conditions on the phase pull that each qubit in- 
duces {(pi = tan^^ (t^)) become stricter if we demand 
a unique determination of the state of the qubits from 
the phase measurement. It is not sufficient that each 
phase pull be different, but rather an additional condi- 
tion is required for unique determination. This is, that 
no phase can equal the sum of any of the others (mod 
2tt). For large N, this could conceivably be a difficult 
constraint, and could make such a cavity-mediated qubit 
measurement scheme demanding. 

In this subsection, we assume that this constraint can 
be met and numerically examine the performance of the 
adaptive scheme described above for an example constel- 
lation with a large number of phases to distinguish. We 
consider four qubits coupled to the cavity mode, and thus 
have 16 different phases to distinguish. The individual 
phase pulls of each qubit are: = 7r/16, 02 = 'I'/S, (/)3 = 
7r/4, (/)4 = 7r/2. The 16-element constellation is shown in 

Fig.m 

Again, we compare the adaptive scheme introduced 
above to static heterodyne detection (here the phase was 
cycled at SOOtt rads/sec in order to be above the thresh- 
old mentioned above) . This comparison is summarized in 
Fig. [51 As in Fig. Oc we show the average a posteriori 
probability (averaged over 500 runs, and for a particu- 
lar choice of correct phase) as a function of the probe 
signal magnitude for three fixed times. As in the small 
N case, the two schemes are almost indistinguishable at 
early times. However, for larger times the performance 
of the adaptive scheme becomes much better than that 
of the static heterodyne scheme. This example suggests 
that the difference in performance between the static and 
adaptive strategies will become more pronounced as the 
number of phases to distinguish becomes larger. The fig- 
ure also shows that the discrimination takes longer as TV 
increases (for both the static and adaptive strategies). 

Like Figure [51 Figure [51 shows the results for a par- 
ticular correct phase. In Table [TTl we show the average 



performance of Vcorr for |a| = 5 and two time values, av- 
eraged over all possible correct values of the phase. Note 
that, as in TableHl there are two levels of averaging here, 
once over many runs of the discrimination task and once 
over the 16 possible correct states. As in the illustrative 
case described by Fig. [51 we see that the average perfor- 
mance of the adaptive scheme is consistently better than 
that of the static scheme. 



Time (sees) 


Detection scheme 


Average Pcorr 


a 


0.2 


Static 


0.3304 


0.0259 


0.2 


Adaptive 


0.4109 


0.0334 


1 


Static 


0.6195 


0.0430 


1 


Adaptive 


0.7994 


0.0442 



TABLE II: = 16 detection: values of the average a poste- 
riori probability of the correct phase, and standard deviation 
of this quantity, at a short and long time instant at fixed SNR 
(|a|=5). 



FIG. 7: The phase constellation (for \a\ = 1) for the discrim- 
ination task described in section IIII B 21 (The dots simply 
show the centers of the coherent states, the uncertainties in 
phase space are not shown). 



IV. CONCLUSION 

We have studied the use of adaptive homodyne mea- 
surements for the task of phase discrimination. This task 
is especially relevant to quantum computing architec- 
tures such as circuit-QED that use harmonic oscillator 
modes to mediate the measurement of qubits. We have 
shown that adaptive measurements provide no perfor- 
mance advantage in the task of measuring one qubit (dis- 
tinguishing between two phases). Following this, we pre- 
sented numerical evidence from comparison of an adap- 
tive homodyne strategy with a static heterodyne detec- 
tion scheme that indicates an advantage in using adap- 
tive homodyne measurements when the task is to mea- 
sure more than one qubit. This advantage is two fold: 
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FIG. 8: (Color online) Contrasting the performance of the 
adaptive and static schemes for discriminating the 16-element 
constellation of section UlI B 21 This figure shows the average 
a posteriori probability of the correct hypothesis as a function 
of probe signal magnitude (SNR) a at three fixed times. The 
dotted (blue) lines plot results for the static scheme and the 
solid (red) lines plot results for the adaptive scheme. As in 
Fig. [5jb) , the error bars at each data point show the standard 
error of the mean values plotted. 

a decrease in the signal-to-noise ratio (which is essen- 
tially the magnitude of the measurement probe signal) 
required to perform the discrimination task to a given 
accuracy, and a decrease in the time required to perform 
the discrimination task to a given accuracy. 

Although we demonstrated the utility of our scheme 
using the circuit-QED setup, it applies equally well to 
other qubit-oscillator systems. In particular, such an 
adaptive scheme will likely be advantageous in improving 
the distinguishability of the states of a single spin when 
measured by coupling the spin to microcavities [13, [H] ■ 

One could imagine extending this study in a number 
of ways. The first would be to consider optimality in 



the N > 2 discrimination task. This could be done by 
formulating a suitable differentiable performance func- 
tion amenable to standard techniques of optimization, 
or by applying more exotic optimization techniques to a 
stochastic performance function. Second, the two exam- 
ples studied in section IIII Bl suggest that the advantage 
offered by an adaptive scheme increases as the number 
of phases to distinguish is increased. It would be inter- 
esting to determine exactly how the advantage (of this 
adaptive scheme over the best static scheme for the con- 
stellation) scales with the constellation size, N. Finally, 
we have not addressed issues of inefficient measurement 
here; all the above calculations have assumed unit effi- 
ciency measurement. Inefficiencies in the measurement 
would modify the POVM of Eq. PH)) . and therefore the 
expression for the likelihood functions and a posteriori 
probabilities. We expect that measurement inefficiency 
will have the same effect as reducing the SNR (probe 
signal size, a), however, this remains to be confirmed. 

Another interesting avenue for further study is the util- 
ity of adaptive schemes for qubit state discrimination 
when the cavity dynamics is non-linear. Such non-linear 
probe devices have recently been developed [10, |3l| , and 
the non-linearity of the cavity results in a more respon- 
sive measurement of the qubit. An investigation of the 
applicability of adaptive discrimination techniques, such 
as the ones described in Refs. [l^l and [13], to such de- 
vices would be interesting. 
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